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1 Introduction 

Several authors have noticed the signature of log-periodic oscillations prior 
to large stock market crashes [[jj, @, |J. Unfortunately good fits of the 
corresponding equation to stock market prices are also observed in quiet 
times. To refine the method several approaches have been suggested: 

• Logarithmic Divergence: Regard the limit where the critical exponent 
(3 converges to 0. |J 

• Universality: Define typical ranges for the free parameters, by observ- 
ing the best fit for historic crashes. Q 

We suggest a new approach. From the observation that the hazard-rate in 
j|] has to be a positive number, we get an inequality among the free variables 
of the equation for stock-market prices. 

Checking 88 years of Dow-Jones-Data for best fits, we find that 25% of those 
that satisfy our inequality, occur less than one year before a crash. We 
compare this with other methods of crash prediction, i.p. the universality 
method of Johansen et al., which followed by a crash only in 9% of the cases. 

Combining the two approaches we obtain a method whose predictions are 
followed by crashes in 54% of the cases. 
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2 The hazard rate 



In pi Johansen et al suggest, that during a speculative bubble the crash 
hazard rate h(t), i.e, the probability per unit time that the crash will happen 
in the next instant if it has not happened yet, can be modeled by by 

h(t) m B (t c - ty a + Bi(t c - t)~ a cos(wlog(t c - t) + ip). 

By assuming that the evolution of the price during a speculative bubble 
satisfies the martingale (no free lunch) condition, they obtain a differential 
equation for the price p(t) whose solution is 



los (f))^{' !(t> ' 



before the crash. Here k denotes the expected size of the crash. 

This implies that the evolution of the logarithm of the price before the crash 
and before the critical date t c is given by: 

(*) log(p(t)) «p c - ^B (t c - tf - K B x {t c - ty 9 cos(o;log(t c -t) + 
P VP + ^ 

With /3 = 1 — a, p c the price at the critical date, and <fi a different phase 
constant. 

Now the hazard rate is a probability and therefore positive. This leads to a 
necessary condition: 

< h(t) 

< B (t c -t)~ a + Bi(t c - t)~ a cos(wlog(t c - t) + ip) 
< B + Bi cos(u log(t c -t) + tp) 

since t < t c . At some times near the critical date cos(w log(t c — t)+ip) takes 
on the values —1 and 1. This implies the necessary conditions 

< B ±B t <^ \Bt\ < Bo- 
On the other hand these conditions are also sufficient for h(t) > since 
cos(lj log(t c — t) + ip) is always between —1 and 1. 

To summarize, if the assumptions of Johansen et al are valid, we must have 
I -Si I < Bq prior to a critical crash. 



3 88 Years of Dow Jones 

To check this model of speculative bubbles, we have investigated the Dow 
Jones index from 1912 to 2000. This period contains 23668 trading days. 
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Johansen et al define a crash as a continuous drawdown (several consecutive 
days of negative index performance) larger than 15%. The following diagram 
shows the drawdowns of the Dow Jones index from 1912 to 2000. Observe 
that there have been 4 drawdowns larger than 15% namely in the years 1929, 
1932, 1933 and 1987. 




As our basic data set, we have calculated numerically the best fit of equation 
(*) to a sliding window of 750 trading days, every 5 trading days. This yields 
4761 best fits. The complete data set is available at |ittp://btm8x5.mat.uni- 
bayreuth.de/~bothmer . 

In what follows, we will call a crash prediction successful, if it was issued 
at most one year before a crash. With this definition there are 229 best fits 
that could possibly give a successful crash prediction. By predicting crashes 
randomly, one would obtain a successful prediction in 4.8% of the cases. 
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3.1 Mean square errors \ 



The first approach to detecting speculative bubbles is to look for good fits 
of (*) to the Dow Jones Index. If the mean square error \ °f the fit is 
sufficiently small one issues a crash prediction. 

The next figure shows the mean square errors of all our best fits compared 
with the best fits before a crash. Unfortunately small errors also occur in 
quiet times. 
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If one issues a crash prediction if the mean square error is smaller then 0.75 
one obtains: 





before crash 


not before crash 


no crash prediction (\ > 0.75) 


175 


2799 


crash prediction {\ < 0.75) 


72 


1732 



I.e. only 72/(72 + 1732) « 3.9% of the predictions are successfull. Since 
this is worse than issuing random predictions, we conclude that one can not 
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predict a crash by looking only at the mean square error. This observation 
has also been made by Sornette and Johansen. 

3.2 Critical Times t c 

If the model of Johansen et al is correct one should expect that a crash occurs 
close to the critical date t c . Using this one can issue a crash prediction when 
the critical date t c is less than one year away. Using our dataset this lead 
to: 





before crash 


not before crash 


no crash prediction (t c > today + 1 year) 


93 


2652 


crash prediction (t c < today + 1 year) 


136 


1879 



I.e 136/(136 + 1879) « 6.7% of the predictions are successful. This is slightly 
better that random predictions, but still not very good. 



3.3 Universality 

Johansen et al suggest that speculative bubbles exhibit universal behavior. 
This would imply that (3 and uj take on roughly the same values for each 
speculative bubble. The following diagram shows the distribution of uj before 
crashes compared with the distribution during other times. 
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One can clearly observe an unexpected peak around u = 9 before the 
crashes. If we issue a crash prediction in the range 7 < uj < 13 we ob- 
tain: 





before crash 


not before crash 


no crash prediction 


150 


3728 


crash prediction (7 < uj < 13) 


79 


803 



I.e 79/(79 + 803) « 8.9% of the predictions are successful. 
The distribution of (3 before crashes and not before crashes is: 
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Distribution of Beta 

Dow Jones 1912-2000 




beta 

Here we observe a tendency toward lower values of j3, but no clear peak. We 
interpret this as evidence, that one should look for logarithmic divergence, 
i.e. the limit of f3 tending to 0, as suggested by Vandewalle et al. ||. We 
will investigate this approach in a later paper. 

3.4 Positive hazard rate 

In section |2] we have explained that in the model of Johansen et al. the 
hazard rate h(t) must be positive. We proved that this is equivalent to 

|#i| < B Q . 

From our best fits we can calculate the value 

k(Bo-|5i|), 

which should also be positive during a speculative bubble, since k is a posi- 
tive number. Consequently we can issue a crash warning if k(Bo — \B\\) is 
positive. With our dataset we obtain: 
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before crash 


not before crash 


no crash prediction (k(Bq — \B\ |) < 0) 


133 


4255 


crash prediction (k(Bq — \B\\) > 0) 


96 


276 



I.e 96/(96 + 276) 25,8% of the predictions are successful. This is already 
a practical success rate, but we can do even better, if we combine this with 
universality. 



3.5 Positive hazard rate and universality 

Combining the last two approaches we issue a crash prediction, if the hazard 
rate is everywhere positive and uj is in the range of section |3.3j . This gives 





before crash 


not before crash 


no crash prediction 


164 


4476 


crash prediction 
(k(B - l-Bil) > and 7 < u < 13) 


65 


55 



I.e 65/(65 + 55) m 54.1% of the crash predictions are successfull. 

The following diagram shows when these crash predictions where issued. For 
every trading day we have plotted the number of crash predictions during 
the past year and the drawdown of the Dow Jones index. 
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Notice that the crashes of 1929 and 1987 have been predicted well in advance. 
The crashes of 1932 and 1933 have not been directly predicted, but we argue 
that they are in the aftermath of 1929 and represent the bursting of the 
same speculative bubble. The crash predictions of 1997 where followed by 
two small crashes in October 1997 and 1998, which didn't quite reach 15%. 
One could argue that they represent a crash in two steps. 



4 Summary 

We have derived a new restriction of the free variables in the model of 
Johansen et al Q for stock market prices during a speculative bubble. This 
restriction alone yields crash predictions with a 25% successrate for the Dow 
Jones index. This is an improvement over the 9% successrate obtained by 
using universality. Combining our approach and the universality method we 
obtain a success rate of 54%. 

We think that these results represent strong evidence for the model of Jo- 
hansen et al describing speculative bubbles in the stock market. 
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